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Abstract 

We test the performance of a recently proposed fundamental measure density functional of aligned hard 
cylinders by calculating the phase diagram of a monodisperse fluid of these particles. We consider all pos- 
sible liquid crystalline symmetries, namely nematic, smectic and columnar, as well as the crystalline phase. 
For this purpose we introduce a Gaussian parameterization of the density profile and use it to minimize 
numerically the functional. We also determine, from the analytic expression for the structure factor of the 
uniform fluid, the bifurcation points from the nematic to the smectic and columnar phases. The equation of 
state, as obtained from functional minimization, is compared to the available Monte Carlo simulation. The 
agreement is is very good, nearly perfect in the description of the inhomogeneous phases. The columnar 
phase is found to be metastable with respect to the smectic or crystal phases, its free energy though being 
very close to that of the stable phases. This result justifies the observation of a window of stability of the 
columnar phase in some simulations, which disappears as the size of the system increases. The only impor- 
tant deviation between theory and simulations shows up in the location of the nematic-smectic transition. 
This is the common drawback of any fundamental measure functional of describing the uniform phase just 
with the accuracy of scaled particle theory. 

PACS numbers: 64.70.Md, 61.20.Gy, 05.20.Jj 
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I. INTRODUCTION 



Monte Carlo simulations conducted on systems of hard anisotropic particles (spherocylinders 
being the most paradigmatic shape) showed that the purely entropic nature of hard core interac- 
tions is enough to explain the stability of different liquid-crystalline phases and phase transitions 
between them lli |2i |3|]. These phases, in decreasing order of symmetry, are known as isotropic 
(I), nematic (N), smectic-A (Sm), columnar (C) and crystal (K) — the isotropic and the crystal 
not being liquid crystalline phases properly speaking—, and some of their physical and chemical 
properties have been described in detail in Refs. Later, Monte Carlo simulations were also 

employed to calculate the full phase diagram of fluids of freely -rotating hard spherocylinders ^] 
and hard-cut spheres [7], including non-uniform phases as the periodic one-dimensional (Sm), 
two-dimensional (C) and three-dimensional (K) phases. 

Several density functional theories have been devised to determine the phase behavior of the 
hard sphere (HS) fluid. These theories can be grouped in two different sets. The first one, the 
weighted-density functionals, are constructed from the knowledge of the thermodynamical and 
structural properties of the uniform fluid [ sl,!^ 1^, while the second one, the fundamental measure 



functionals (FMF), initially introducec 
description of the HS freezing I isi 



by Rosenfeld Jll , Il2|] and later improved for an adequate 



1511 . are built on the geometry of the particles alone. 



The extensions of these theories to hard anisotropic particles have not been as successful as 
they have been for HS. There are two reasons to explain this difficulty: The first one is related to 
the, as of today, still poor knowledge of the structural properties of fluids composed by anisotropic 
particles, and the second one is the inherent complexity in dealing with orientational degrees of 
freedom within density functional theory. This notwithstanding^ome weighted-density function- 
als have been developed for the fluid of hard spherocylinders ivl to study both the I-Sm and 
the N-Sm phase transitions as a function of the particle aspect ratio. These functionals were con- 
structed as modifications of a reference HS weighted-density functionals, and their predictions, 
tested against Monte Carlo simulations, are reasonably good. They do not allow though to prop- 
erly account for the C and K phases. 

FMF are more appropriate to treat these phases as, by construction, they are more suitable 
to describe highly confined particles, such as they are in a solid. Unfortunately the fundamental 
measure formalism has little flexibility to apply it to arbitrary geometries. FMF have been obtained 



for parallelepipeds with restricted orientations of their principal axes [|18l 



and very recently 



2 



for cylinders also with a parallel alignment constraint piQ. For freely rotating anisotropic 
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particles 



25|],but 



FMF have been obtained for needles, infinitely thin plates, and their mixtures [|22| , 
this time the price to pay is to eliminate at least one of the characteristic lengths of the particles. 
Besides, the numerical minimization of these functionals to obtain the equilibrium density profiles 
of non-uniform phases seems to be a very demanding task. 



J] 



In this article we aim at testing the recently proposed FMF for parallel hard cylinders 
by comparing its predictions with Monte Carlo simulations reported in the literature [1261 
We will consider all possible non-uniform phases, namely N, Sm, C and K and will depict the 
phase diagram the FMF predicts. There is an interesting aspect about this model that poses a 
particularly stringent test on the theory. In Ref. t2M a window of stability of the C phase was 
reported whose existence the authors of Ref. 1270 could not completely settle, although their results 
pointed to its being a finite size effect because this window disappears — being preempted by a 
K — in simulations of very large systems. We will show that our FMF does indeed confirm this 
conclusion by showing that either the Sm or the K are always more stable than the C, although the 
difference in free energy is rather small — what justifies its observation in small systems. We will 
also compare the resulting equations of state for the N, Sm and K phases with those obtained from 
the Monte Carlo simulations of Ref. [12711 and conclude that the performance of our functional is 
almost perfect in the description of highly non-uniform phases, even improving on the free- volume 
description of the K phase. 



II. FUNDAMENTAL MEASURE DENSITY FUNCTIONAL 



In [12 ih we obtained a fundamental-measure density functional for mixtures of parallel hard 
cylinders, so we will just gather here the formulae, specialized for the case of a one-component 
fluid. The functional is constructed out of the one for two-dimensional hard disk. There are 



two versions for the latter: Rosenfeld's original version il2|l . and the version of Tarazona and 
Rosenfeld nl5\\ . The former has some important drawbacks, for instance, the low density limit of 
the functional is only approximate. That of Tarazona and Rosenfeld recovers the exact result in this 
limit. On the other hand, the former is easier to implement than the latter, because it is expressible 
in terms of one-particle- weighted densities, while that of Tarazona and Rosenfeld contains a two- 
particle-weighted density. Nevertheless both are amenable to numerical treatment and we will 
explore the results of both. So the formulae presented here will describe the implementation of the 
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two versions for the functional of parallel hard cylinders. 

Irrespective of the version we are using, the free-energy density functional can always be writ- 
ten 

/?^[p]=/?J-id[p]+/?^ex[p], (1) 

where (3 is the inverse temperature in units of the Boltzmann constant, 

/?^id[p] = JdrJ dzp{r,z) [lnVp(r,^) - 1] (2) 

is the functional of the ideal gas (V is the thermal volume, irrelevant for the phase behavior), and 
PJ^cx[p] is the excess free energy due to interactions. We are using the notation r = {x,y) for 
vectors perpendicular to the cylinders axes. Fundamental-measure functionals are expressed in 
terms of an excess free-energy density <l>(r, z), such that 

/3^e.[p] = jdrj dz<l>{r,z). (3) 

This free-energy density can be given as a function of a set of weighted densities. The whole set 
of them can be written in terms of the two densities 

Po(r, z) = ^ [p(r, z + L/2) + p(r, z - L/2)] , (4) 

rz+L/2 

pi{r,z) = / p{r,t)dt. (5) 

Jz-L/2 

Common to both versions are the weighted densities 

1 

2^ 



^o{r,z) = / po{r + K,z)dR, (6) 



ni{r,z) = I pi(r + R, z) dR, (7) 



|R|=i? 
|R|=i? 

n2ir,z) = J po{v + R,z)dR, (8) 

|R|<_R 

ns{r,z) = J pi{r + R,z)dK. (9) 

|R|<_R 

For Rosenfeld's original version [12] there are also two vector densities, namely 

vi(r,2;) = ^ j po(r + R,2;)RciR, (10) 

|R|=iJ 

V2(r,2;) = ^ j pi(r + R,^)RdR, (11) 

|R|=iJ 



riiria + 27ri?^ nprii - vi ■ V2 2 ^1 - vj 
$Ros = -'^o ln(l - ng) H ^ h vri? ^2 7- r^. 



and the expression for the excess free-energy density is 

(12) 

For Tarazona-Rosenfeld's version [|15|] there are also two two-particle-weighted densities, namely 



|Ri|=iJi IR2NR2 
1 - 

'2R 



N2{t,z) = jdRi J dR2Piir + Ri,z)piir + R2, z)K (^^^^^^y (14) 



|Ri|=iii IR2NR2 

where 



K{x) = — Vl — x'^ sm~^ X, (15) 



and the expression for the excess free-energy density is 



$ = -noln(l -na) + — + - (16) 



III. PHASE BEHAVIOR 

The Euler-Lagrange equation 

provides the equilibrium density for the system when there is no external field and the chemical 
potential is fixed to p (equivalently, when the mean density is fixed to the value p corresponding 
to that chemical potential). Expected phases are nematic (no spatial ordering), smectic (one- 
dimensional layering of particles), columnar (two-dimensional odering of liquid columns) and 
crystal (a combination of both orderings). These are the phases shown in the simulations of Veer- 
man and Frenkel [27]. Quite as expected, columnar phase is a triangular ordering of columns and 
crystal phase is a piling up of such triangular lattices, i.e. what is commonly referred to as an AAA 
crystal (see Fig.[T])- 

A direct solution to (flTI) is numerically unfeasible so, as it is customary, we have resorted to 
a variational method. Thus, in order to account for all the above phases in our density functional 
description in a unified simple way, we have chosen the parametrization 

P(r) = pKTx±ir)xiiiz), (18) 



FIG. 1 : (Color online) Triangular (AAA) crystal. The lattice parameters a and d are shown in the figure. 



where p is the mean density (number of particles per unit volume) and 

X±(r) = ^ ^exp [-a^ (r - K^f] , (19) 

k 

/aw \ 1/2 

The parameter v}^i is defined as the D-dimensional volume of the unit cell of the corresponding 
phase (D = 1 smectic, D = 2 columnar, D = 3 crystal). Its values are 

= d, KiiJ = v^«V2, = dV3a'/2, (21) 

d being the layer spacing along the Z direction and a the lattice parameter of the triangular lattice 
on the XY plane (see Fig.[T]). Finally, Rk = kiai + k2a2 (ki, k2 G Z), with a„ = ^ (^Vs, (— l)"j 
the vectors defining the two-dimensional triangular lattice. In Appendix |X] we give explicit ex- 
pressions for the weighted densities evaluated with the density profile (fTSi) . 

When Eq. (flTl) . using the parametrization (fTSi) . leads to a solution with a\\ = a_L = 0, the 
equilibrium phase is a nematic; a smectic is the equilibrium phase if ay ^ and a± = 0; it is a 
columnar if ay = and a± ^ 0; and a crystal if both a\\ ^ and a.^ ^ 0. For the crystal phase 
1 — pV^cu = y provides the fraction of vacancies. 
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A. Nematic phase 



When a\\ = = in (fTSl) . both (fT2l) and (fT6l) provide the same free-energy density, namely 

_ BFv — 

^ = ^ = % + r]{\ny + 3y + y''), (22) 

where $o = ?7ln(V/f ) — r/, a linear term irrelevant for phase behavior, -q = pv h the packing 
fraction, v = -kR^L is the volume of a cylinder, and y = r]/{l — r]). This free-energy density is 
plotted in Fig.|2l 

From (|22|) the equation of state is readily obtained as 

I3pv = y + + 2y^ = r] (23) 



the same equation of state as that of parallel hard cubes [|28l] . 

The structure factor can also be obtained from the relationship S{q, q^) = [1 — pc(g, ^z)]"^, 
where c(g, g^) is the Fourier transform of the direct correlation function of the uniform fluid. Its 



expression was given in Ref. \2M [cf- Eqs. (39)-(43) and Appendix B]. Specializing to the one- 



component fluid and taking into account that 



j d(\5{u 


— r) = 


2TTU'^o{qu), 


(24) 


j dqe{u 


— r) = 


iru'^^i^qu), 


(25) 


J dq,Qiu/2- 


-kl) = 


u^2iqzu/2), 


(26) 


|q|, r = |r|, *o(a;) 


= M^) 


, ^i{x) = 2Ji{x)/x and \l'2(a;) 


= sinx/x. 



Jo{x) and Ji (x) being the zeroth and first order Bessel functions, respectively, we obtain, from the 
Tarazona-Rosenfeld functional (fT6l) . the following expression for the inverse structure factor 



=1 + Sy^,{2q*)m2{2ql) + V [2^ oiq*)"^ M*)"^ 2{2q:) + vl/i(2g*)vl/2(g:)'] 

+ 2y2(l + 2y) [2vI/o(g*)*i(g*)^2(g:)' + vl/i(g*)2vl/2(2g:)] (27) 
+ 1/2(1 + 62/ + 62/2)i[/,(g*)2x[/2(g:)2, 

where g* = Rq and g* = Lqz/2. 



1 



_Q 4 I 1 1 1 1 1 1 1 1 1 

'0.3 0.4 0.5 0.6 0.7 

FIG. 2: Free-energy densities = $ — $0 — air/— a2 (with oi = 4.8463 and 02 = —2.0555 chosen so as to 
amplify the differences between the different free-energy branches) vs. packing fraction rj for the nematic 
(solid line), smectic (dashed line), columnar (dotted line) and crystal (dash-dotted line) phases. The N- 
Sm bifurcation point is shown by a filled square. The nematic-columnar and smectic-columnar coexisting 
packing fractions are marked with open squares and open circles, respectively. The columnar phase is 
metastable and hence so are these two phase transitions. Before the smectic changes into a columnar the 
crystal becomes more stable. The smectic-crystal phase transition is marked with full circles. 

B. Smectic phase 

When we set = in (fTSl) and substitute this density profile into either (fT2)) or (fT6l) . both 
yield the same expression 

m = n.{z) |- In [1 - n,{z)] + + "^^f | , (28) 

i - n3[z) [1 -n3{z)\ ) 

with 

rioiz) = ^ [p{z - L/2) + p{z + L/2)] , (29) 

rz+L/2 

n^lz) = ttR^ / dz'p{z'). (30) 

Jz-L/2 
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So both theories predict the same nematic-smectic transition. 

Solving Eq. (fTT]) . a solution with a\\ 7^ is found for every r] > r/N-sm ~ 0.31 (also plotted in 
Fig.|2l). The value of a\\ approaches zero as i] approaches this value from above. On the other hand, 
the free-energy density for this smectic phase is tangent to that of the nematic one (see Fig. O, 
so the transition is continuous. This being so, we can obtain a more accurate value of r/N-sm as 
the smallest 77 at which the structure factor (|27]) diverges at some wave vector q = 0, 7^ 0. 
Specializing (|27]) for such a wave vector we find 

S{0, q,)-' = 1 + 2y{A + 5y + 2y^)^2{2q:) + y\9 + Uy + 6i/2)^2(g:)'. (31) 

The smallest r] for which the right-hand side of (|3TI) vanishes at a g* is r^N-Sm = 0.3143, and the 
value of g* at which it happens corresponds to a smectic period d/L = 7r/g* = 1.3015. 

C. Columnar phase 

At packing fraction ?7^_c — 0.4369 the nematic loses stability against columnar ordering. This 
value is determined from the divergence of the structure factor (|27l) at a wave vector q 7^ 0, ^2 = 0, 
which, for the Tarazona-Rosenfeld functional (fT6l) . is given by 

S{q,0)-' = l + 4y{2 + y)^i{2q*)+4y\3 + 2y)^o{q*)'^i{q*)+y\S + 10y + Qy')^l{q*). (32) 

In this case, however, the columnar free energy is not tangent to the nematic one, so the transition 
is first order. We can determine the N-C coexistence by the usual double tangent construction. This 
yields the ?7n = 0.3957 and r]c = 0.4425 as the coexisting packing fractions of the nematic and the 
columnar phases, respectively (see Fig.[2l). At the latter, the lattice parameter is a/R = 2.4744. 

We can see here an important difference between this version of the functional and that based 
on Rosenfeld's original approximation, Eq. (fT2)) . The latter leads to the following equation for the 
inverse structure factor 

S{q,0)-' = l + 2y(2 + y)(vI/^(g*)-|*o(g*)r)+2y(2 + 7y + 42/2)^o(g*)^i(g*) 

+y'{3 + Wy + 6y')^l{q*), (33) 

where the new complex vector ^o{q*) = iJi{q*)c{/q has been introduced. The value of ?7j^_c 
which this approximation predicts is r^^-c ~ 0.5599. If we had to believe this value for the N-C 
bifurcation, the columnar free energy would be much too high to be consistent with the metastable 
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columnar phase found in simulations 1261 2711 . For this reason, we have not pursued this version 
of the funcional anymore. 

The columnar free energy is higher than the smectic one up to = 0.6534, where a first 
order Sm-C transition occurs, with coexisting packing fractions ?7sm = 0.6382 and r]c = 0.6697. 
However at these packing fractions the equilibrium phase is no more the smectic but the crystal, 
thus the columnar phase is always metastable, and in particular so are the N-C and the Sm-C 
transitions. All this can be easily visualized in Fig. |2l 



D. Crystal phase 

At packing fractions around r] ^ 0.58 a solution to Eq. (fTTl) with «|| ^ and a± ^ renders 
a free energy smaller than that of the, up to that point stable, smectic phase. The fluid undergoes 
a first order Sm-K transition with coexisting packing fractions r/sm = 0.5689 and r/K = 0.5936. 
The lattice parameters of the coexisting crystal are a/R = 2.3102 and d/L = 1.1419. With these 
values the fraction of vacancies can be found to be just a mere 0.3%. The crystal is the only stable 
phase for rj > rj^up to close packing (see Fig.©. 



IV. COMPARISON WITH COMPUTER SIMULATIONS 



Numerical simulations for this fluid were carried out first by Stroobants et al. [26] and later by 
Veerman and Frenkel [|27ll . The former, made with 900 cylinders, showed the sequence of stable 
phases N-Sm-C-K. The latter confirmed this result but also made simulations with 1080 cylinders 
which showed that the columnar phase previously found appeared due to a finite size effect. Their 
conclusion was that the columnar phase is always metastable, but has a free energy very close to 
that of the smectic phase, so much that the boundary conditions may artificially render it more 
stable. Our previous calculations are fully consistent with this result, as Fig. |2]illustrates. 

Besides this first qualitative agreement, we can also perform a more quantitative comparison 
with simulations by comparing the equations of state. This is done in Fig.[3l The simulation results 
are those obtained with the largest system size 1270. The figure shows that the agreement between 
the numerical values of the pressure is excellent for all stable phases. The values for the crystal 
phase are indistinguishable from the simulations, as it is also the location of the Sm-K transition. 

The only important deviation between theory and simulations concerns the location of the N- 
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Sm transition. While both, theory and simulation, predict that this transition is continuous, the 
theory predicts that it occurs ati] = 0.3143 while the simulations yield a value ofr] = 0.443. This 
failure of the theory to predict the location of continuous transitions between low-density uniform 
and non-uniform phases is a fingerprint of FMT. For instance, the FMF of parallel hard cubes 
also predicts the same value of rj = 0.3143 for the transition between the fluid and the smectic, 
columnar and crystal phases (the later being the stable one) HQ, while simulations provide a 



value of ?7 = 0.49 for the freezing of this fluid IS; 
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3lh . The reason for this drawback lies in the fact 



that, by construction, FMFs provide, in the uniform limit, the SPT equation of state — which for 
anisotropic bodies deviates from the exact result — , while at the same time the prediction for the 
nonuniform phases improves significantly due to the dimensional crossover properties of FMFs 



||32D- This discrepancy in the accuracy with which the theory describes both type of phases leads 
to inaccurate predictions of the uniform- nonuniform phase transition points. 

We end this section by comparing the EOS for the crystal phase given by the FMF and that 
obtained by a cell approximation for the fluid of parallel hard cylinders, which is derived in Ap- 
pendix |Bl Figure |4] shows the results of both theories as well as the simulation results. As it can 
be seen, while the FMF results fit perfectly the simulation points, the cell approximation, although 
still a rather good description, underestimates the EOS. We can also see that, as expected, both the- 



ories converge at high densities, a 



crossover 3D-^0D of the FMF UJ, 



mown result which is a direct consequence of the dimensional 



V. CONCLUSIONS 

There are very few examples in the literature in which the same functional describes with ac- 
curacy all inhomogeneous phases of a liquid crystalline fluid. In this article we have applied a 



fundamental-measure functional recently proposed for mixtures of parallel hard cylinders [|21|] to 
determine the phase behavior of the one-component fluid. As usual with fundamental-measure- 
based functional, the results obtained for the uniform (nematic) fluid are those provided by scaled 
particle theory, and so the accuracy the functional provides for this phase is reasonably good but 
not perfect. As a consequence, the predicted nematic-smectic phase transition deviates signifi- 
cantly from the Monte Carlo simulations of Refs. Iz^H]. although the order is correct. However, 
the accuracy with which the remaining stable phases, smectic and crystal, are obtained is excellent, 
the plots being indistinguishable from the simulation data, even for the smectic-crystal coexisting 
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> 

CQ. 




FIG. 3: Equations of state (reduced pressure vs. packing fraction) for all stable phases obtained from the 
fundamental measure functional for parallel hard cylinders. These phases are: nematic (for packing frac- 
tions up to the point indicated by a full rhombus), smectic (from that point up to the discontinuity) and 
crystal (from the discontinuity up to close packing). The open circles are the simulation results reported in 
Ref. [27]. Arrows mark the nematic-smectic and smectic-cystal phase transitions as obtained from those 
simulations. The two insets show the equations of state for the columnar metas table phase in the neigh- 
borhood of the nematic-columnar (left inset) and smectic-columnar (metastable) phase transitions. [Labels 
stand for nematic (N), smectic (Sm), columnar (C) and crystal (K).] 

densities. Results for the equation of state of the crystal improve on those obtained by a cell 
approximation (which we have also reported in an appendix). Another correct prediction of the 
theory is that the columnar is only a metastable phase, but its free energy is sufficiently close to 
that of the stable phases so as to justify the observation of a window of stability of that phase in 



the oldest simulations 1261 ] made with the smallest system size, a window that disappears when 



the size in increased [27]. In summary, the proposed functional provides excellent results, very 
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FIG. 4: Comparison between the equation of state of the crystal phase as obtained from minimization of the 



functional (solid line), from the cell approximation (dashed line) and from simulations 12711 . 



similar to those obtained by simulations, but obtained at a much cheaper price. They also made us 
confident that its version for mixture may provide very good results as well. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR THE WEIGHTED DENSITIES 

Insertion of the parametrization (fTSl) into the expressions for the weighted densities dil)-® 
leads to the formulae 



no{r,z)=pV^SQ^1\v)pf\zl 
n,{r,z)=pV^^S,^Q^'l\r)Hf\z), 
n2{r,z)=pV}^,^T^f\r)pf\z), 
ns{r,z) = pV}£,^T^f\r)Hf\z), 



(Al) 
(A2) 
(A3) 
(A4) 



where V^^^'' is defined in Eq. (|2T]) . The functions are given in terms of 

1/2 _2 , , 1 



9a{x) 



ea{x) = - erf (a/ox), 



erf (x) being the standard error function. To be precise, 



(A5) 



Q^'Iir) = 1, 

Q^'lir) = Q(5(r) = (7aJi?)5^^?aJ|r-Rk|)/o(2i?ax|r-Rk|), 



(A6) 
(A7) 



where In stands for the zeroth-order modified Bessel function of the first kind. The rest of the 



(A8) 

T^Hr) = 27r V(7,^(|r-Rk|) / dt t g^^{t) Io{2ta^\r - R^l), (A9) 
k -^0 

(AlO) 

I i^-kd + L/2) + g^^^ {z - kd - L/2)], (Al 1) 



expressions are similar; 




= 7rR\ 




= T^l\r) 




= 1, 




= Pfiz) 


Hf{z) 


= L, 




= Hf{z] 



Hf (z) = J2 ban {z-kd + L/2) - e„|| {z-kd-L/2)]. 



(A12) 
(A13) 



As for the two-particle weighted densities, after a lengthy calculation (see Ref. [|2l|l for some 
details) Ni can be expressed as 



N,{r,z) = 2{pV}S,^fpf\z)H^i'\z)J^^\r 



(A14) 
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with the functions P^^^^ and H^^^ defined above. The radial contribution is 

/j^(r) = ttR^, (A 15) 

^ ki,k2 

PIT 

X / dttsint Io[Bi^-^^i^^it,r)], (A16) 
Jo 

where 

5ki,k2(t, r) = 2i?a^^Cki + Ck2 + 2CkiCk2 cos(t + - ^kj, (A17) 
denoting r - Rk, = Ck,(cos^/'k„, sin V'kJ, with u = 1,2. Finally, N2, is given by 

N2{r,z) = [pV}^,^Hf\z)r/P{r). (A18) 

APPENDIX B: CELL APPROXIMATION FOR THE CRYSTAL PHASE OF PARALLEL HARD 
CYLINDERS 

This section is devoted to obtain a cell approximation for the free energy per particle of the 
crystal phase of parallel hard cylinders. To this aim we first calculate the free volume available to 
one particle moving in an cell defined by the first nearest neighbours: a prism with hexagonal base 
composed by six triangular cells of period a (see a sketch in Figure [5]) and height equal to 2d. Six 
hard disks (the cylinder sections) of radii R are fixed at the vertexes of the hexagon while a seventh 
one is allowed to move within this cell, with the only constraint of not overlapping the other six 
disks (which of course do not overlap themselves). Simple geometric considerations lead, for the 
area accessible to the center of mass of the seventh disk, to the formula 



V3x^ + cos ^ X — x\J\ — 

3 . 



(Bl) 



where x = a/AR. The free volume of this cell is simply Vf^ee = '^^freeL{y — 1) with y = d/L. 
If we fix the mean packing fraction of the crystal, the variables x and y are related through the 
equation rj = v/V^^-^^, where v = txR^L and V^^l is defined in (|2TI) . are the particle and cell 
volumes respectively. Thus y = n/sVSrjx'^. 

The cell theory approximates the free energy per particles as 

^ = -ln(^), (B2) 
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FIG. 5: (Color online) Sketch of the tiiangular lattice of period a. The free region of area Afrcc within which 
one particle can move appears colored. One sixth of this area can be obtained substracting from the area of 
the triangle ABC those of the triangle ABD and of the sectors BED and ADF. 



with V the thermal volume of the system, which in our case is 



In ( ^ I - In f VSx^ + cos"^ x - xVl - - - In ( — J 1 ) . (B3) 



Once the mean packing fraction is fixed the free-energy (IB3I) must be minimized with respect to x 
with the constraint x > 1/2 (x = 1/2 represents the close packed limit), and then the pressure is 

obtained as (3Pv = r^^T^, with the result 

or] 



l-Axlr]/r]^p 



Vcp = tt/v 12 being the value of rj at close packing, and Xq the solution to the equation 

(cos-^ x-^^+x (^VSx - Vl - x2 j = 0. (B5) 
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